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ABSTRACT 

Context. Near-infrared (hereafter NIR) data may provide complementary information to the traditional optical population synthesis 
analysis of unresolved stellar populations because the spectral energy distribution of the galaxies in the 1-2.5 /im range is dominated by 
dilferent types of stars than at optical wavelengths. Furthermore, NIR data are subjected to less absorption and hence could constrain 
the stellar populations in dust-obscured galaxies. 

Aims. We want to develop observational constraints on the stellar populations of unresolved stellar systems in the NIR. 
Methods. To achieve this goal we need a benchmark sample of NIR spectra of "simple" early-type galaxies, to be used for testing 
and calibrating the outputs of population synthesis models. We obtained low-resolution (R~1000) long-slit spectra between 1.5 and 
2.4/jm for 14 nearby early-type galaxies using SofI at the ESO 3.5-m New Technology Telescope and higher resolution (R~3000) 
long-slit spectra, centered at the Mgl at ~ 1 .5 1 fim for a heterogeneous sample of 5 nearby galaxies observed with ISAAC at Antu, one 
of the 8.2-m ESO Very Large Telescope. 

Results. We defined spectral indices for CO, Nal, Cal and Mgl features and measured the strengths of these features in the sample 
galaxies. We defined a new global NIR metallicity index, suitable for abundance measurements in low-resolution spectra. Finally, we 
present an average NIR spectrum of an early-type galaxy, built from a homogenized subset of our sample. 

Conclusions. The NIR spectra of the sample galaxies show great similarity and the strength of some features does correlate with the 
iron abundance [Fe/H] and optical metal features of the galaxies. The data suggest that the NIR metal features, in combination with 
a hydrogen absorption feature may be able to break the age-metallicity degeneracy just like the Mg and Fe features in the optical 
wavelength range. 

Key words. Infrared: galaxies — Galaxies: abundances — Galaxies: elliptical and lenticular, cD — Galaxies: spiral — Galaxies: 
stellar content 



?H 1. Introduction often lead to non-unique solutions for the stellar populations. 
^ „ „ , , , • ■ , , In this respec t , the Lick /IDS system of i ndices pioneered by 
State-of-the-art space-based mstrumentation can only resolve iBursteiri e t al.'(^198?) and'Faber et al.'flQSS), and developed fur- 
galaxies front the Local Group well and in the more distant ob- ^ ^^^^^^ , ^^^g^ 2OOO, 2005), has been particularly 
jects we usually can see only the tip of the red giaiit branch stars ^^^^^^^f^^ for interpreting the integrated optical light of galax- 
which IS rarely surhcient tor population analysts. This leaves us 

withthedifficiilttaskoftryingtorecoverthestellarpopulationof However, new constraints are necessary to interpret more 
unresolved galaxies froin their integrated properties. The com- Seated systems and one possibility is to widen the spectral 
plex mix of stellar populations found in most of them usually ^^^^.^^ ^j^^ near-infrared (NIR) because the light in dif- 
makes It possible to constrain only the most recent generation of ^^^^^^ wavebands is dominated by different populations of stars, 
stars. However, such properties as the stellar kinematics and the ^j^^ passbands are dominated by light from older and red- 
present-day metal content are dominated by the overall star for- therefore offer us the possibility to study other stellar 
ination histo ry, and toget h er with the well known age-met^ populations than is possible with optical spectra alone. In ad- 
degeneracy mSmi O'Connelll liM ^fflS they ^^^^ determinations through optical spectroscopy 
Send ojfprint requests to: M. Cesetti P^^^i^l^ heavily reddened evolved stellar populations 
* Based on observations made with ESO Telescopes at the La Silla ^^'^^ as dusty spheroids or some bulge globular clusters hidden 
and Paranal Observatory under programmes 69.B-0311 and 077.B- by dust (e.g., [ Stephens & Frogel| 2004). NIR spectroscopic ob- 
Qlg3 servations could overcome these problem because the extinction 
** Table 5 is only available in electronic form at the CDS i" the /T-band is only one-tenth of that in the V-band. 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via Most of the previous work at NIR was focused on ei- 
[http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-i- A/ 1 ther active galactic nuclei (AGNs) or objects with very strong 
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Star formation, including recent surveys of ultra-luminous in- 
frared galaxies ('Goldader et^[1995t iMurphy et al." 1999', '200 1^, 
Burston et al. 2001), lu minous infrared galaxies (Goldader et al. 
1997; Reunanen et alj l2007h . starbursts (Engelbracht 1997; 
Coziol et al. 200V), Seyfert galaxies (Ivanov et al. ,2000: 
Sosa-Brito et al. 2001; Reunanen et al. 2002; Boiss onet al.l 
2002), LINERS (Larki n et al.ii 19981: lAlonso-He rrero et al."2000'; 
Sosa- Brito et al. 2001), and interacting galaxies ( Vanz i & Rieke 
19971: IVanzi et al.lll998h . 



Relatively few NIR spectroscopic observa tions exist for 
"normal" galaxies. Only iMannucci et al.l (1200 ll) provided low- 
resolution (R ~ 400) template spectra for galaxies of different 
Hubble types, including some giant ellipticals. Such data, to- 
gether with the corresponding analysis is a necessary first step 
towards developing a system of spectral diagnostics in the NIR 
because well-understood galaxies with relatively simple star 
forming history will all ow us to tune the NIR population syn- 
thesis models. Recently. [Silva et alj (l2008l hereafter SOS) stud- 
ied the stellar populations of eleven early-type galaxies in the 
nearby Fornax cluster by means of A'-band spectroscopy. 

A few promin ent NIR features were first studied by 
lOriglia et al.l (Il993h who demonstrated that they represent a su- 
perb set of indicators for constraining the average spectral type 
and luminosity clas s of cool evolved s t ars. Th is co nclusion was 
later c onfirmed by iForster Schreibej ( |2000|) and llvanov et al.l 
(12004 . Furthermore, the s ame NIR fea tures appear to be promis- 
ing ab undance indicators dFrogel etani200ll: [Stephens & Frogell 
l2004l) . 

Here we describe two new data sets of high quality NIR spec- 
tra of ellipticals/spirals designed to provide a benchmark for fu- 
ture NIR studies of unresolved galaxies: (i) low -resolution spec- 
tra covering the range from 1 .5 to 2.4 fim that include strong fea- 
tures such as CO, Nal and Cal that are traditionally studied, and 
( ii) moderate resolution spectra around the Mgl absorption fea- 
ture at 1 .5 1 fim. This is the second strongest Mg feature in the H- 
and /T-band atmospheric windows (after the feature at 1.71 fim) 
and at zero redshift it is located in a region of poor atmospheric 
transmission, making it difficult to observe in stars. However, the 
redsh ift of externa l galaxies moves it into a more transparent re- 
gion (Ilvanovl l2001), just the opposite of the 1.71 jjm Mgl which 
becomes affected by the red edge of the //-band atmospheric 
window. 

We are only few years from the launch of the James Webb 
Space Telescope (Gardner et al. 2006), a space-based infrared 
telescope that will have unprecedented capabilities. Therefore, 
in the near future the application of NIR spectroscopy to the 
study of galaxy properties will be limited not by lack of data but 
by our understanding of spectral features at these wavelengths. 
Improving the characterization of the NIR indices is a timely 
step in this direction. 

The paper is organized as follows. The sample selection is 
discussed in Sect.|2] The NIR spectroscopic observations and 
data reduction are described in Sect.[3] The definition of the new 
NIR spectral indices and their measurements are given in Sect.H] 
Results are discussed in Sect. |5] and summarized in Sect.|6l 



ellipticals (NGC4478, NGC4564, NGC4742, NGC5077, 
NGC 5576), one is a bright dwarf elliptical (NGC 3641), and two 
are bulges of disk galaxies (NGC 4281, NGC 4594; our spec- 
tra are dominated by their spheroi ds). The Mg2 index spans the 
range between 0.21 and 0.34 mag dBenderet al.ll 1993b . in agree- 
ment with other Mg2 literature values, and the central velocity 
dispersions are between about 130 and 370km s"^ The galaxy 
properties are summarized in Table[T] 

Four galaxies in the sample exhibit evidence for weak nu- 
clear activit y: NGC 4- 5 94 an d NGC 5077 were identified as 
LINERs bv iHeckmanI (Il980l) ; NGC 4472 is a weak Seyfert 2 
(Ho et al. 1997) and IC4296 shows radio and X-ray emission 
( Rinn et al. 2005) . Finally, NGC 4649 is a close companion of 
the giant spiral NGC 4647 a nd it may have undergone some tidal 
stripping (lDas& Joglll999l) . 

The low resolution used for the observation of our sample is 
not sufficient for characterizing the 1.51 yum Mgl feature, so we 
selected from the ESO Science Archive higher resolution spectra 
covering the region around ffie 1.51 fim for a heterogeneous set 
of galaxies. Although these observations were obtained for other 
purposes, ffiey provide us with a first glimpse into the behavior of 
this feature. The properties of these galaxies are also described 
in Table[l] 



2. The sample 

Our main sample consists of 14 nearby bright (Bj- <14.1 mag) 
and undisturbed spheroids. It was selected to cover a wide 
range in luminosity and velocity dispersion and to have available 
Lick/IDS line-strength indices, from the literature. Six of our 
galaxies are giant ellipticals (NGC 4472, NGC 462 1 , NGC 4649, 
NGC 4697, NGC 6909, IC4296), five are intermediate-size 



Table 1. Basic properties of the sample galaxies. 



ID Galaxy 


R.A.(2000) 


Decl.(2000) 


Type 


Redshift 


cr 


Bj 


H.J 


Kj 


Age 


Mg2 


Mg2 




Fe5335 


[Z/H] 


[Fe/H] 


CaT 












(kms-') 


(mag) 


(mag) 


(mag) 


(Gyr) 


(mag) 


(mag) 


(A) 


(A) 


(A) 


(A) 


(A) 


(1) (2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12a) 


(12b) 


(13) 


(14) 


(15) 


(16) 


(17) 












Galaxies with low-resolution spectra 
















1 IC 4296 


13 36 38.85 


-33 57 59.3 


E 


0.0124(a) 


inn/ \ 

337.0(a) 


11.61 


14.56 


14.30 


5.2(j) 


0.323 


0.340(a) 


1.22(m) 


3.03(m) 


0.39G) 






2 NGC3641 


11 21 08.85 


+03 1 1 40.2 


Epec 


0.0058(b) 


163.2(i) 


14.10 


15.92 


15.17 


6.6(k) 


0.284 


0.273(k) 


0.72(k) 


2.39(k) 






6.343 


3 NGC3818 


11 41 57.50 


-06 09 20.0 


E5 


0.0055(a) 


187.5(a) 


12.67 


14.69 


14.54 


6.4(k) 


0.315 


0.322(a) 


1.67(k) 


2.87(k) 


0.37(o) 


0.16(o) 


6.242 


4 NGC4281 


12 20 21.52 


+05 23 12.4 


SO+:sp 


0.0089(c) 


230.5(h) 


12.25 


13.91 


13.80 




0.314 














5 NGC4472 


12 29 46.76 


+07 59 59.9 


E2 


0.0033(a) 


309.7(a) 


9.37 


15.34 


15.01 


lO.O(k) 


0.306 


0.331(a) 


1.37(k) 


3.11(k) 


0.26(o) 


0.06(o) 


6.119 


6 NGC4478 


12 30 17.53 


+ 12 19 40.3 


E2 


0.0047(a) 


153.0(a) 


12.36 


15.62 


15.45 


6.9(k) 


0.253 


0.260(a) 


1.67(k) 


2.81(k) 


0.30(o) 


0.16(o) 


6.810 


7 NGC4564 


12 36 27.01 


+ 11 26 18.8 


E 


0.0038(a) 


171.1(a) 


12.05 


14.74 


14.66 


8.3(k) 


0.321 


0.329(a) 


1.54(k) 


2.81(k) 






6.325 


8 NGC4594 


12 39 59.43 


-11 37 22.9 


Sa(s) sp 


0.0034(a) 


259.1(a) 


8.98 


12.22 


12.22 


9.5(k) 


0.330 


0.338(a) 


1.44(k) 






0.39(o) 


6.361 


9 NGC4621 


12 42 02.39 


+ 11 38 45.1 


E5 


0.0014(a) 


232.0(a) 


10.57 


10.28 


10.16 


9.1(k) 


0.328 


0.345(a) 


1.39(k) 


2.98(k) 


0.07(o) 




6.055 


10 NGC4649 


12 43 40.19 


+ 11 33 08.9 


E2 


0.0036(a) 


368.5(a) 


9.81 


14.61 


14.35 


11.9(m) 


0.338 


0.347(a) 


1 29(k) 


2.57(k) 


0.29(o) 


0.06(o) 




11 NGC4697 


12 48 35.70 


-05 48 03.0 


E6 


0.0041(a) 


171.7(a) 


10.14 


13.35 


13.31 


6.4(k) 


0.297 


0.296(a) 


1.69(k) 


1 94fk) 


0.06(o) 


-0.03(o) 


6.185 


12 NGC5077 


13 19 31.66 


-12 39 25.0 


E3-4 


0.0090(d) 


255.9(i) 


12.38 


16.21 


14.93 


15.0a) 


0.295 


0.324(j) 


1.83(m) 


2.64(m) 


0.12(i) 






13 NGC5576 


14 21 03.60 


+03 16 14.4 


E3 


0.0050(e) 


183.1(i) 


11.85 


14.94 


14.74 




0.253 


0.243(f) 


1.52(n) 


3 11(1) 




0.02(o) 




14 NGC6909 


20 27 38.60 


-47 01 34.0 


E+: 


0.0090(f) 


128.7(f) 


12.61 


14.25 


14.32 




0.208 


0.214(f) 


2.17(p) 


2.32(p) 


















Galaxies with high-resolution spectra 
















15 Mrkl055 


02 48 18.55 


-08 57 37.4 


S? 


0.0366(v) 






14.69 


14.30 


















16 NGC 1144 


02 55 12.32 


-00 11 01.7 


E 


0.0288(q) 




13.78 


15.18 


15.06 


















17 NGC 1362 


03 33 53.08 


-20 16 57.7 


SO": pec 


0.004 l(t) 


91(s) 




16.11 


15.68 








0.223(f) 










18 NGC 4472 
















see above 


















19 NGC 7714 


23 36 14.28 


+02 09 17.8 


Sb(s):pec 0.0093(p) 


62(u) 


13.00 


15.85 


16.08 
















4.7(u) 



Columns: (1) Identification number, used in the plots, (2) Galaxy name, (3-4) Coordinates, (5) Morphological type from lde Vaucouleurs et al. I ( 1 19911 . hereafter RC3), except for Mrk 1055 and 
NGC 1 144 f rom the Lyon Extragalactic Database (LEDA), (6) Redshift, (7) Central velocity dispersion, (8) Total observed blue magnitude from RC3, (9)-( 10) Total obs e rved H and K magnitude 
from lCutri e t al. (20031, 2 MASS All-Sky Catalog of point sources), (11) Age, (12)-(17) Optical indices available in literature, (12a) the Mg2 values from lBender et alj ( |1993|) and (17) Calcium 
Triplet from Cenarro et al. (2003). 

Refere nces: (a)- Smith et al. (2OO0|), (b)- lFalco et all jl999h . (c )-lBinggeli et al.l jl985l) . (d)- lda Costa et al.l jl998|) . (e)-lDenicol6 et al.l j2005h . (f) -IWegner et"aLl (l2()03h. (g)- lKobavashi & Arimotd 
dl999h. (h)J alle Ore et al. (1991), (i)-Faber et al. (1999), (j)-"Annibali et al. (2007), (k)- Sanchez-Blazauez et al. (2007), (!)■ Howell (2005), (m)-RamDazzo et al.' ('2005'), (n)-Trager et al. (1998), 
(o)- lTrager et al, C200a) . (p)- .Ogando et al.. C2008) , (q)- .Keel (1996) , (r) -.Lauberts & Valentiin (.1989) , (s) -.Bernardi et ah (.2002) . (t)- ,01iva et al.. a995.) . (u)- .Gardner et al.. t2006) . (v)- ,Lu et al., C1993.) . 
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3. Observation and data reduction 

3.1. Soft low-resolution spectra 

Long-slit (0'.'6 x 290") NIR spectroscopy of 14 galaxies was 
perfor med with the Sofl (So n of ISAAC) infrared spectrom- 
eter (Mo orwood et alJ[T998bl) at the Nasmyth A focus of the 
ESO 3.5-m New Technology Telescope (NTT) at the European 
Southern Observatory (ESO) in la Silla (Chile). The detector 
is 1024x1024 pixels Hawaii HgCdTe array with 18.5-/zm pixel 
size. It is read out in four quadrants and the average quantum ef- 
ficiency is ~65%. We used the low resolution grism (R ~ 1000) 
providing spectral dispersion of 10.22 A pixel"' and a coverage 
between 1.53 /vm and 2.52 fim. 

The smaller galaxies were observed in nodding mode, plac- 
ing them at two different positions along the slit, separated by 
150". In the case of bigger galaxies, the telescope pointing was 
alternated between the galaxy and a nearby clear sky region to 
obtain a good sky sampling. We randomly dithered the telescope 
within a 10" box along the slit around each position. This dit her- 
ing, also known as jitter (see e.g. Sofl handbook or iDevillardI 
11999.) . helps us to remove bad pixels, to improve the pixel sam- 
phng, and the flat-field correction because on each image the 
object is placed on slightly different position on the array. 

Immediately before or after each galaxy spectrum we ob- 
served a B star to measure the atmospheric transmission and the 
instrument response, always in nodding mode and with a simi- 
lar airmass to that of the corresponding galaxy. Finally, to cali- 
brate the spatial distortions we also obtained a sequence of stellar 
spectra placed on 19 equidistant positions along the slit. 

The details about the spectroscopic observations are given in 
Table|2] 

The spectra were reduced with the Image Reduction and 
Analysis Facility (IRAFfl The sky emission was removed either 
by subtracting the corresponding image in each nodding pair or 
by subtracting from each object image the average of the pre- 
ceding and the succeeding sky images, respectively for the two 
modes of observation described above. The data were flat-fielded 
with screen flat spectra. Geometric distortion correction was ap- 
plied along both the spatial and dispersion axes by using an arc 
frame and the grid of stellar spectra mentioned above, respec- 
tively. Simultaneously, the frames were wavelength calibrated. 
We applied the dispersion corrections, and removed the ho t pix- 
els an d cosmic ray hits with the algorithm described in iPvchI 
(|2004|) . Then, we aligned the individual 2-dimensional spectra 
along the spatial direction (i.e., along the slit), to combine them 
together and to extract a 1 -dimensional spectrum of the galaxy 
with apertures as listed in Table|2] 

Finally, we removed the telluric absorption. To do this we 
first created a sensitivity function by multiplying the observed 
spectrum of the standard by the intrinsic spectrum of a star with 
the corr esponding (or the closest available) spectral type from 
iPicklesI (il998) . Some residual emissions were observed, indi- 
cating that the spectral types of the standards may not be ac- 
curate. We fitted these artificial features with a Gaussian, with 
the IRAF task SPLOT and subtracted them from the sensitiv- 
ity function. This prompted us to omit the bluest part of the 
spectra (A <1.65jum) from further analysis because the stan- 
dards are rich of Hydrogen lines that modify the appearance 
of the galaxy spectra. Next, we divided the galaxy spectra by 



Table 2. Log of the spectroscopic observations 



ID 


Galaxy 


Obs. Date Exp. 


Slit Galaxy 


Standard Star 


Aper- 


S/N 






YYYY- Time 


P.A. 


secz 


Sp.Type 


sec z 


ture 








MM-DD (sec) 


(deg) 








(arcsec) 




(1) 


(2) 


(3) (4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 






Galaxies with low-resolution spectra 








1 


IC4296 


2006-04-17 1200 


0.0 


1.19 


B2V 


1.19 


6.8 


85.5 


2 


NGC3641 


2006-04-17 3600 


70.0 


1.20 


B2V 


1.21 


4.0 


60.2 


3 


NGC3818 


2006-04-17 3600 


100.0 


1.34 


B2IV 




5.7 


81.2 


4 


NGC4281 


2006-04-17 2400 


178.0 


1.26 


B2V 


1.26 


4.0 


73.5 


5 


NGC4472 


2006-04-16 2000 


155.0 


1.32 


B5 


1.32 


17.0 


83.7 


6 


NGC4478 


2006-04-16 2400 


140.0 


1.39 


B5V 


1.39 


8.5 


64.5 


7 


NGC4564 


2006-04-16 1800 


137.0 


1.57 


B5 


1.57 


5.7 


89.2 


8 


NGC4594 


2006-04-16 1200 


0.0 


1.42 


B5V 


1.42 


5.7 


54.7 


9 


NGC4621 


2006-04-17 1200 


165.0 


1.33 


B2IV 


1.33 


5.7 


85.2 


10 


NGC4649 


2006-04-16 2000 


105.0 


2.27 


B5III 


2.27 


22.7 


88.7 


1 1 


NGC46y7 


2006-04-16 1400 


70.0 


1.33 


B5 


1.33 


8.5 


85.7 


12 


NGC 5077 


2006-04-18 3200 


10.0 


1.24 


B8IV 


1.24 


6.8 


75.2 


13 


NGC 5576 


2006-04-16 1500 


95.0 


1.65 


B2.5V 


1.65 


5.7 


102.7 


14 


NGC 6909 


2006-04-17 2000 


70.0 


1.09 


B2V 


1.09 


2.8 


52.5 






Galaxies with high-resolution spectra 








15 


Mrk 1055 


2002-08-07 1800 


179.5 


1.05 


G3V 


1.04 


2.9 


78.0 


16 


NGC 1144 


2002-08-09 1800 


117.7 


1.21 


G3V 


1.68 


1.5 


52.0 


17 


NGC 1362 


2002-08-19 1800 


179.5 


1.02 


B3V 


1.23 


1.5 


67.0 


18 


NGC 4472 


2002-05-19 900 


179.5 


1.22 


B3III 


1.54 


5.9 


71.0 


19 


NGC 7714 


2002-07-19 1800 


179.5 


1.14 


B3III 


1.03 


2.3 


96.0 



(1) Identification number, used in the plots, (2) Galaxy name, (3) 
Observing Date, (4) Total exposition time, (5) Slit position angle, (6) 
Airmass, (7) Standard star spectral type, (8) Standard star airmass, (9) 
Length of the apertures used to extract 1 -dimensional spectra shown in 
Figs |2|3| andll II (10) Average signal-to-noise per 1/pixel (see ®. 



this sensitivity function to restore the overall continuum shape 
of the galaxy spectrum. This process, broken into steps, is il- 
lu strated in Fig. [ T| and mor e detailed descr iptions can be found 
in lHanson etaU (Il996l) and lMaiohno et alj (.1996) . The final re- 
duced spectra were taken to zero redshift (the values used are 
listed in Table[lJ and are shown in Figs.|2]and[3] 



' Image Reduction and Analysis Facility is a general purpose soft- 
ware system for the reduction and analysis of astronomical data, IRAF 
is written and supported by the IRAF programming group at the 
National Optical Astronomy Observatories (NOAO). 



3.2. ISAAC high-resolution spectra 

Five spectra of galaxi es were obtained with IS AAC - a NIR 
imaging spectrometer jMoorwood et alJll998al) at Nasmyth B 
focus of Antu, one of the 8.2-m ESO VLT unit telescopes of ESO 
Very Large Telescope at Paranal (Chile). It is equipped with a 
1024x1024 Hawaii Rockwell HgCdTe array with a spatial scale 
of0.146"pixeLi. The long-slit (a.'3 x 120") medium (R~3000) 
resolution spectroscopy mode was used. The wavelength cover- 
age spans the range from ~1.48jum to ~1.62/vm and the central 
wavelength was chosen according to the object redshift. A sum- 
mary of the observations is given in Table|2] 

Similar to the Sofl observations, telescope pointings ware 
alternated between the galaxy and a nearby empty sky region, 
or if the dimension of the galaxy were small enough, then by 
nodding among two positions along the slit separated by about 
20". No jittering was used. Standards were observed in nodding 
mode soon before or after the galaxy, at a similar airmass. A 
sequence of stellar spectra placed at equidistant positions along 
the slit was obtained to calibrate the spacial distortions, as in the 
case of the Sofl data. 

The reduction steps were carried out using IRAF and they 
are the same as for the main data set, albeit with some complica- 
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NGC 4478 extracted spectrum 
before telluric correction 



Observed spectrum of 
Ktandard star (Hip 024825} 




(b) 



Reference spectrum of the same NGC4478 extracted spectrum 
spectral type of the standard star (D3 V) after telluric correction 



X 2 





MOOO 20000 

WavGlcngth{A ) 



26000 14000 20000 26000 

¥avclcngth(A ) 



Fig. 1. The telluric absorption correction: the wavelength cal- 
ibrated 1 -dimensional galaxy spectrum (a) is divided by the 
wavelength calibrated 1 -dimensional observed spectrum of the 
standard star (b) and the result is multiplied by the reference of 
th e same spectra l type of the standard star (c) from the library 
of lPicklesI (Il998h . The product (d) is the corrected galaxy spec- 
trum. Note the large scatter at ~18000-19000/im region where 
the Earth's atmosphere is practically opaque to the radiation. 



tion due to the lower signal to noise and longer exposure times. 
First, the individual images were cleaned from cosmic rays using 
SPEC_COSM, a cosmic ray re jection script for long slit spec- 
tra developed by Ivan DokkumI d2001,) . Second, bad pixels were 
masked out before the combination because the small number of 
images limited the ability of the rejection algorithms to remove 
the bad pixels during the combination of the individual images. 
Finally, the highly variable sky background led to residuals in 
the sky-subtracted images which were removed with the IRAF 
task BACKGROUND before the geometric distortion correction 
and the wavelength calibration. 



4. Definition of the NIR spectral indices 

The NIR spectral indices are usually defined with a central band- 
pass covering the spectral feature of interest, and two other 
bandpasses at the red and blue sides which are used to trace 
a local continuum le vel through a lin e ar fit to the m ean values 
in both bands, as i n IaU et all (11995), Frogeletal. (2001) and 
llvanov et al.l (l2004l) . When this is not possible, for example be- 
cause the feature is at the end of the atmospheric window, then 
either one continuum band or a fit to the overall shape of the 
continuum is used, as in Kl einmann &l!ail ([19861) . 

Most of the earlier definitions are aimed at stars, where the 
spectral features are typically only a few kms"' wide, however 
in galaxies the intrinsic velocity dispersion broadens them to a 
few hundred kms This forced us to define new spectral in- 
dices with broader passbands (Table[3] Fig.|2]and[3]l. The con- 
tinuum passbands were placed on regions clear of absorption 
features, and not (or only weakly) affected by the telluric absorp- 
tion. Furthermore, we tried to have as wide continuum passbands 
as possible, to increase the signal-to-noise of the measurements. 
The passband of each spectral feature was defined to encompass 
the maximum extent of the feature itself and to avoid, if possi- 
ble, contamination by other features. For instance, the passband 



i 2 - 




l.SOxlO* 1.55x10* 1.60x10* 1.65x10* 1.70x10* 1.75x10* 
A (A ) 

Fig. 2. //-band spectra of the sample galaxies, continuum- 
divided (i.e., all continua were normalized to 1) and shifted verti- 
cally for displaying purposes by adding increasing shifts of 0.25. 
The position of the main spectral features is shown with dotted 
lines. The passbands of the Mgl (1.71 fim) feature and its adja- 
cent continuum regions are marked with inclined solid and dot- 
ted lines, respectively. 



of our CO (2.30 fim) index is wide but still only includes the 
'^CO feature and there is no contribution from the '^CO band 
head. 

We measured the equivalent width (EW) of the lines with 
respect to a local continuum obtained by fitting a straight line to 
the "clear" parts of the spectrum: 



EW = (1 -Ffeature//^co„t)*A^ 



(1) 



where /^feature is the flux inside the feature window, Fcont is the 
value of the local continuum linearly interpolated (or extrapo- 
lated, if necessary) at the center of the line and AA is the width 
of the feature bandpass. The fluxes /^feature and Fcom are normal- 
ized by the bandpass widths. The measurements were performed 
following the above recipe with a self-written IDL script. The 
error-bar are derived in the same script by means of Monte Carlo 
simulations on the measured spectrum considering as noise the 
RMS of the normalized flux in the continuum region. The av- 
erage signal-to-noise calculated on the continuum passbands re- 
gions is written in Table |2] The corresponding values of each 
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Nal Cal CO '=CO"COCO 



' IC4296 




NGC4472 



_ NGC4478 



_ NGC45g4 y 
' NGC4621 




_ NGC5077 
' NGC5576 
_ NGC6909 




Table 3. NIR Spectral indices. 

The central wavelength and bandwidth A/1 are given for the cen- 
tral feature (Line) and the continuum bandpass at shorter (Cont. 1) and 
longer (Cont. 2) wavelengths. 



Feature 


Line 




Cont. 1 


Cont. 2 




/I, 




X, 


A.1 


K 


A/i 




(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


Mgl(1.50yum) 


15040 


80 


14960 


80 


15120 


80 


Mgl(1.71 yum) 


17098 


90 


17033 


46 


17206 


100 


Nal (2.21 pim) 


22070 


140 


21965 


62 


22170 


40 


Cal (2.26 yum) 


22647 


106 


22553 


74 


22802 


172 


CO (2.30 /zm) 


22990 


160 


22802 


172 







Table 4. Equivalent widths of the newly defined NIR spectral 
indices for the sample galaxies. 



Galaxy 



Mgl 
EW ±o-(EW) 
(A) 



Nal 
EW±o-(EW) 
(A) 



Cal 
EW±o-(EW) 
(A) 



CO 
EW±D-(EW) 
(A) 



2.1x10' 



2.2x10* 



2.3x10' 
A (A) 



2.4x10* 



IC 4296 
NGC 3641 
NGC3818 
NGC 4281 
NGC 4472 
NGC 4478 
NGC 4564 
NGC 4594 
NGC 4621 
NGC 4649 
NGC 4697 
NGC 5077 
NGC 5576 
NGC 6909 



2.44: 
2.54: 

3.23 : 

3.07 : 
3.04: 
3.50: 
2.97 : 
3.24: 
2.73 : 
2.77 : 
3.14: 
3.72 : 
2.37 : 
3.60: 



0.40 
0.43 
0.20 
0.26 
0.44 
0.29 
0.36 
0.39 
0.23 
0.40 
0.41 
0.56 
0.30 
0.48 



5.74: 

4.48 : 

5.64: 
3.35: 
6.11 : 
4.45 : 
6.69: 
7.35 : 
6.25: 
6.49: 
5.85 : 
1.90: 
4.27 : 
2.59: 



0.37 
0.40 
0.19 
0.28 
0.18 
0.43 
0.15 
0.35 
0.11 
0.16 
0.20 
0.15 
0.15 
0.46 



4.01 ± 0.35 
1.88 ±0.56 
2.67 ± 0.53 
3.08 ± 0.59 
2.72 ±0.17 
1.91 ±0.81 
2.25 ± 0.41 
3.01 ± 0.40 
2.79 ± 0.28 
3.10 ±0.25 
3.10 ±0.28 
2.78 ± 0.84 
1.75 ±0.41 
2.72 ± 0.77 



15.48 
14.28 
13.63 
16.00 
14.01 
14.29 
14.84 
14.65 
15.09 
14.19 
14.88 
14.64 
12.63 
15.22 



±0.63 
± 0.57 
± 0.65 
± 0.80 
± 0.63 
± 1.00 
±0.71 
± 1.12 
±0.64 
± 0.60 
± 0.58 
± 1.06 
± 0.70 
±0.94 



Fig. 3. /T-band spectra of the sample galaxies, continuum- 
divided (i.e., all continua were normalized to 1) and shifted verti- 
cally for displaying purposes by adding increasing shifts of 0.25. 
The position of the main spectral features is shown with dotted 
lines. The passbands of the Nal (2.21 /zm), Cal (2.26 /im), and 
CO (2.30 yum) features and their adjacent continuum regions are 
marked with inclined solid and dotted hnes, respectively. 

equivalent width with their relative errors for all galaxies in the 
sample are given in TableU] 

Our definitions have intentionally wide bandpasses to make 
them intrinsically insensitive to the varying velocity broadening 
of spectral features in the galaxy spectra. To verify this, we mea- 
sured the indices for broadened spectra of M and K type giants. 
We choose stars in this spectral range because they are the most 
representative of the early-type integrated galaxy spectra. (These 
spectra are a reasonable match to the observed galaxy spectra.) 
We computed for them the correction factor C(cr) for various 
line-of-sight velocity distribution, as in SOS: 



C(cr) = 



EW(o- = 0) 
EW(cr) 



(2) 



the EW of the feature is measured from the stars convolved with 
a Gaussian with cr given in brackets. We broadened the stellar 
spectra to velocity dispersion up to 400 km s ' in steps of 50 km 
s"' and used templates for different spectral stars and from dif- 
ferent sources. The corrections for the different spectral indices 
are shown in Fig.|4] Some of them show a significant scatter. 



depending on the adopted stellar spectrum. However, the mean 
correction are typically small (i.e., ^ 15% in the worst case). 



5. Discussion and Conclusions 

5. 7. General appearance of the spectra 

The spectra of our sample galaxies appear qualitatively simi- 
lar in most of the NIR features except for a large spread in 
the Nal values is present. This similarity is not surprising be- 
cause we selected mostly giant ellipticals and spheroids, which 
are all relatively metal rich, have no significant recent star 
formation, and onl y weak nuclear a ctivity , if any. This con- 
clusion agrees with iMannucci et al.l (1200 ll) who demonstrated 
that galaxies within the same Hubble type have nearly identi- 
cal NIR spectra. On the other hand, the large spread of Nal, 
with respect to the observational errors, appears to be real, 
and suggests variety of star formation and enrichment histo- 
ries among the galaxies in our sample. Two of the three galax- 
ies with systematically weaker Nal exhibit strong H/3, indi- 
cating that the weak Nal might be related to the presence of 
younger stellar populations (the third galaxy lacks optical spec- 
troscopy). The strongest features are the CO absorption bands 
in both the H and A'-band atmospheric windows. They origi- 
nate in K and M stars, as can be seen from the libraries of stel- 
lar spectra describe d in Lancon & Rocca-Volmerange (19921), 
IOrigliaetaI.1 (Il993h . IPalher et al., (.1996,) . .Mever et al., (.19981) . 
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Fig. 4. Correction factors of the NIR spectral index for different 
values of velocity dispersion. They have been derived by com- 
paring the EW measured in the spectra of the labeled stars before 
and after broadening them to a particular velocity dispersion. 
The spectral types of the stars are in a range between Kill and 
Mill with a spread of [Fe/H] between -0.35 and -1-0.41. In each 
panel the solid line is the 2nd order polynomial fitting the data. 



[Wallace & Hinkld ( Il997h and lForster Schreibeil (l2000l) . A num- 
ber of weaker metal absorption features are visible as well: Sil 
at 1.589//m, Mgl at l.Tlljum, Nal at 2.206 and 2.209 yum and 
Cal at 2.261, 2.263 , and 2.2 66 /im. They are also present mostly 
in cool stars (i.e. Kleinmann & Hall 1986; Origliaetal. 1993; 
Wallace & Hinkld il997: .Forster Schreiber .2000; Mever et al.. 



19981) . As mentioned above, we only consider the redder fea- 



tures with A>1 .65yum. 

Not surprisingly, our galaxies show no Fell at 1.644yum or 
H2 lines line emission, usually associated with supernova activ- 
ity and supernova related shocks, respectively. There is also no 
Bry emission. 

5.2. Behavior of the NIR spectral indices 

Evolutionary synthesis models (i.e. IWorthevll 1 99 4^ are the usual 
method to interpret the behavior of spectral features in galax- 
ies, but the simple comparison of these features with the galaxy 
parameters can be informative too. For example, the optical in- 
dices, Mg2 , <Fe> and H/j s how tight correlation with the cen- 
tral velocity dispersi on (e.g. Terlevich et al. 198ll: lBernardi et al.l 
Il998t iMehlertet all 12003; Morelli et al. 2008) . suggesting that 
the chemical and the dynamical evolution of ellipticals are inter- 
twined. 

The EW of our NIR spectral indices, the iron abundance and 
metallicity are plotted versus the central velocity dispersion in 
Fig.E] The loci of the S08 data (see their Fig. 13) are shown 
with dashed lines. Among the galaxies, NGC5077 (#12) and 
NGC 6909 (#14) possess very weak Nal, Cal lines, and their rel- 
ative errors are significant. The correlations between the NIR 
indices and velocity dispersion for the sample galaxies are plot- 
ted in FiglH the fit do not consider NGC 4281 (#4), NGC 5077 
(#12) and NGC 6909 (#14) because of the reasons explained fur- 
ther on. The correlation with CO shows a similar slope as the one 
found by S08 but with a different zero-point. The different defi- 
nition of the index and especially the position of the continuum 
bandpass can lead to a systematic variation of the EW values. 
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Fig. 5. The equivalent widths of the NIR spectral indices, iron 
abundance [Fe/H], and metal abundance [Z/H] in the sample 
galaxies as a function of the central velocity dispersion. The la- 
bels correspond to the row numbers in the Table [T] The open 
and filled circles refer to quiescent and active galaxies, respec- 
tively. The small and large circles refer to galaxies with an age 
between 5 and lOGyr and >10Gyr, respectively. Galaxies with 
no known ages were assumed to be old. The dashed lines corre- 
sponds to the relations find by S08. The solid line in each panel 
represents the linear regression (y = ax -i- b) through all the data 
points except for #4, #12 and #14. The Pearson correlation coef- 
ficient (r) and the results of the linear fit are given. Typical error 
bars are shown. 



This effect seems to be more evident in the CO band, whose 
equivalent width is measured with only extrapolated continuum 
blue-wards of the feature. The agreement with S08 is closer in 
the Nal versus cr relation although our data show bigger scatter 
Without pretending to trace a general conclusion we found rela- 
tively large Nal values compared to stars and something similar 
was also found by S08. The iron abundance [Fe/H] and the to- 
tal abundance [Z/H] only show very weak correlations. Active 
galaxies are clustered at the high end of the velocity dispersion 
distribution which is understandable because their velocity dis- 
persion measurement may be affected by the black hole and the 
active nucleus. Note that our highest velocity dispersion galaxy 
is a Type 2 Se yfert and the CO b and may suffer some dilution, 
as discussed in llvanov et al.l (12000 ). 

The EW of our NIR spectral indices are plotted against each 
other in Fig.|6] and compared to the relat ions for cluster star s 
and solar neighborhood giants by S08 and Davi dge et alj ( l2008l) . 
respectively. Three galaxies - NGC 4281 (#4), NGC 5077 (#12) 
and NGC 6909 (#14) - are above the relation of S08 for the clus- 
ter stars, which hints at different chemical enrichment history 
with respect to the rest of the sample. Unfortunately the literature 
lacks much data about these objects and although NGC 6909 has 
the lowest velocity dispersion in our sample, these galaxies do 
not stand out in any respect, including in the optical Mg2 versus 
cr diagram. Further investigation of these galaxies is necessary. 
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Fig. 6. The equivalent widths of the NIR indices plotted against 
each other. The symbols are the same as in Fig.|5] The dashed 
lines correspond to the relation s found by SOS . The dotted line 
corresponds to the relation by Davidge et alJ (,2008) for giant 
stars in the solar neighborhoods. All these relations are plotted in 
their observed range. The solid line represents the linear regres- 
sion (y = ax -H b) through all the data points except for #4, #12 
and #14. The Pearson correlation coefficient (r) and the results 
of the linear fit are given. Typical error bars are shown. 

The rest of our objects populates a locus that follows similar 
trend as the galaxies of SOS. In the Cal versus CO plot we see 
a correlation with significant scatter, and with a different slope 
than in SOS. Finally, the Mgl versus CO plot is dominated by 
scatter. Galaxies with traces of nuclear activity, evident in other 
wavelength ranges than the NIR, do not seem to separate from 
the rest of the sample, and they show no traces of emission lines. 
Therefore, the contribution of their AGNs is negligible with re- 
spect to the rest of the galaxy. 

Figure|7] shows the EW of the NIR spectral indices versus 
the Mg2 measurements from lBender et al.l (Il993h . Nal, and to a 
lesser extent Cal and CO, do show trends with respect to Mg2, 
with NGC42S1 (#4), NGC5077 (#12) and NGC6909 (#14) 
standing out. The correlation of Nal with both cr and Mg2 sug- 
gests that the Nal feature is dominated by the stellar photosphere 
rather than by the interstellar medium. 

Given that [Fe/H] measurements are not available for the en- 
tire sample and that Mg2 is not representative of the total chem- 
ical abundance of a galaxy because of the varying abundance 
ratios, we attempted to create a combined Iron- and-g-element 
index similar to that defined by iGonzalej (119931) and more re- 
cently bv lThomas eTall (l2003h : 

[MgFe]' = ^JMgb X (0.72 x Fe5270 + 0.2S x Fe5335) (3) 

Such a combined indicator is expected to decrease the effect 
from the varying a/Fe ratio. Since not all the components of this 
indicator are available we directly substitute the iron and mag- 
nesium indices by defining a new indicator: 



Fig. 7. The equivalent widths of the NIR indices p lotted against 
the Mg2 indices measured bv lBender et alj ( Il993l) . The symbols 
are the same as in Fig.|5] The solid line represents the linear re- 
gression (y = ax + b) through all the data points except for #4, 
#12 and #14. The Pearson correlation coefficient (r) and the re- 
sults of the linear fit are given. Typical error bars are shown. 



The results are shown in Fig.[8^-c. The Nal index of NGC 5077 
(#12) and the CO index of NGC 6909 (#14) deviate from the 
main loci of the other galaxies. 

To address this issue we plotted the optical indices of these 
galaxies versus the combined [MgFel' index (Fig. [8tl). This is an 
analog of the typical plot (i.e . I Worthevll 1 994) that allows to dis- 
entangle the age-metallicity degeneracy: the inverse H/j index is 
roughly proportional to age while the [MgFe]' is dominated by 
metal abundance. The two galaxies exhibit strong which sug- 
gests that they are dominated by populations of 3 Gyr or younger 
(see for example Fig. 1 in SOS). NGC 5077 and NGC 6909 are 
also well separated from the rest of the galaxies on the versus 
Nal plot (Fig.|8^) which leads us to the conclusion that the NIR 
indices can be used to create a similar diagnostic plot, but mea- 
suring the Bry feature requires better quality data than the ones 
described here. 

We investigated the behavior of the NIR indices in relation to the 
H and K-band magnitude and the H-K color but no clear corre- 
lations were found. 

5.3. Combined NIR metal index 

Ground based NIR spectroscopy is much more time consuming 
than the corresponding optical observations - higher and vari- 
able background and detectors with worse cosmetics often re- 
quire one to sacrifice either resolution or signal-to-noise to ob- 
tain the data in a reasonable amount of time. To alleviate this 
problem, we defined a combined spectral index of all the major 
A'-band metal features: 



[MgFe]" = VMg2 X Fe5335 



(4) 



(CONaCa) = (CO -H Na i -H Ca i)/3 .0 



(5) 
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Fig. 8. (a-d) The equivalent widths of the Nal, Cal, CO and B/i 
indices plotted against the [MgFe]" index, (e) The equivalent 
widths of the H/3 index as a function of the Nal index. The sym- 
bols are the same as in Fig.|5] The dashed lines correspond to 
the relations found by SOS. The solid line represents the linear 
regression (y = ax + b) through all the data points except for 
#4, #12 and #14. The Pearson correlation coefficient (r) and the 
results of the linear fit are given. Typical error bars are shown. 



Various weighting schemes were tried to minimize the scatter of 
the basic relations (Fig.|9]l. However, the simple average yielded 
the tightest relations. The carbon and oxygen are a elements, 
while the sodium and calcium originate in both high and low 
mass stars. The average value of the CO index for the sample 
galaxies is ~ 15 A, the Nal is ~4.3A and the Cal is ~2.1 A which 
means that at least 2/3 of the new index is dominated by metals 
produced mostly in high mass stars, i.e. early in the history of 
the elliptical galaxies. The other implication is that a relatively 
limited amount of recent star formation could affect the new in- 
dex more than it would an iron peak dominated index. The lack 
of iro n peak features in the NIR is well known as it was pointed 
out in llvanovl ( l2001h . 

The new index improves the correlations. For example, 
the Pearson correlation coefficient is ~10% higher for the 
(CONaCa) vs. Mg2 and even -20% higher for the (CONaCa) vs. 
[MgFe]", with respect to the respective correlations where only 
one IR index is used. We derived these relations using only the 
galaxies with known parameters, i.e. excluding the poorly stud- 
ied NGC4281(#4) and, NGC5077(#12) and NGC6909(#14) 
which systematically show peculiarities with respect to the bulk 
of our sample. 

5.4. Average NIR spectrum of the sample galaxies 

Studies of composite stellar systems (i.e. galaxies hosting an 
AGN) often need to subtract the contribution of the underlying 
galaxy. This prompted us to create an average spectrum of the 
galaxies in our sample. We used a homogenized subset of eight 
galaxies, excluding the objects with young populations. We also 



Fig. 9. The equivalent widths of the newly defined <CONaCa> 
as a function of the central velocity dispersion (a), Mg2 (b), H/3 
(c), and [FeMg]" (d) indices. The symbols are the same as in 
Fig-E] The solid line represents the linear regression (y = ax H- 
b) through all the data points except for #4, #12 and #14. The 
Pearson correlation coefficient (r) and the results of the linear fit 
are given. Typical error bars are shown. 

Table 5. Average spectrum of the sample galaxies observed at 
NTT in arbitrary units and normalized to unity at 22000 A. The 
regions with zero rm.s. at the beginning and the at the and are 
covered only by one or two spectra. The full table is given only 
in the electronic version of the journal. 





F, 


r.m.s. 


14944 


2.217 


0.000 


14945 


2.384 


0.000 


14946 


2.384 


0.000 


14947 


2.386 


0.000 


14948 


2.130 


0.366 


14949 


2.572 


0.254 


14950 


2.575 


0.251 


14951 


2.578 


0.248 


14952 


2.583 


0.245 


14953 


2.588 


0.243 


14954 


2.595 


0.241 


14955 


2.602 


0.240 



excluded those with low signal-to-noise ratio. The average age 
of the galaxies used to create the composite spectrum is 9+2 Gyr 
(the median is 9.5 Gyr) and the average [Fe/H] is 0. 17±0. 13 (the 
median is 0. 16). The composite spectrum is shown in Fig.fTOland 
the values of the spectrum are listed in Table|5] together with the 
rm.s. values per Angstrom. We included the //-band section be- 
cause the artifacts caused by the spectral type mismatch of the 
standards are minimized by the averaging of galaxies observed 
at different redshifts. 



10 



Cesetti, M. et al.: Early-type galaxies in the near-infrared: 1.5-2.4 spectroscopy 



c 

D 

h 2.5 



3 




a - 



1.5x10" 



1.6x10" 1.7x10" 
X(A) 



o o o 




2x10" 2.1x10" 2.2x10" 2.3x10" 2.4x10" 
X(A) 



2.5 



E 1.5 






Mrkl055 



NGC1144 



362 



jA^-A^-iiyw^ NGCl 



^A/iAV'^Mll^r^^^^^V^^ NGC4472/M49 



NGC7714/Mrk538 



1.49x10" 1.5x10" 1.51x10" 1.52x10" 1.53x10" 
X (A) 



Fig. 10. Average H- (a) and K-hand (b) spectrum of sample 
galaxies normalized to unity at 22000 A (thick line). The thin 
lines correspond to the + 1 cr confidence levels. Some of the more 
prominent spectral features are marked. 



5.5. The Mgl feature at 1.51 



llvanov et al.l (|2004|) pointed out the possibility of using the Mgl 
feature at 1.51 //m as an a-element abundance indicator in the 
NIR. With exploratory purposes we obtained spectra of a small 
and diverse set of galaxies (Table|2] and Fig.fTTTi. To carry out 
quantitative analysis we defined an index (TableO and for the 
five galaxies we measured equivalent widths of 3.3, 5.0, 3.8, 
3.9 and4.5A, respectively for Mrk 1055, NGC 1 144, NGC 1362, 
NGC4472 and NGC 77 14. The typical uncertainty is ~ 0.3A. 
No corrections for velocity dispersion were applied. Given the 
diverse nature of the objects, we refrain from drawing any con- 
clusions but we note that the relation between the optical and 
the NIR Mg features is not straightforward because the only two 
Mg2 values that we have from the literature are inversely pro- 
portional to our measurements for NGC 1362 and NGC 4472. 
Nevertheless, these observations prove that it is feasible to mea- 
sure the NIR Mgl feature at 1.51 yum and give a basis for future 
NIR synthetic spectral modeling. 

6. Summary 

1. We report on a new data set of NIR spectra of ellipti- 
cals/spheroids. The stronger features were measured and 
show mild correlations with the optical metal features. NIR 
versus NIR features show a lot of scatter but the diagram of 
CO versus Nal seems to be a good diagnostic for detecting 
different chemical enrichment histories. 

2. The NIR metal absorptions, in combination with a Hydrogen 
absorption feature may be able to break the age-metallicity 
degeneracy, but further investigation and better data are 
needed to verify this possibility. 

3. A new combined NIR index is defined that allows one to in- 
terpret lower signal-to-noise data. It is dominated, at about 



Fig. 11. Spectra of the galaxies observed with ISAAC in the re- 
gion of the Mgl (1.51 jum). The central wavelengths of the indi- 
vidual Mgl lines are shown with vertical solid lines. The pass- 
bands for the feature and the adjacent continua are shaded. The 
spectra were normaUzed to unity and shifted vertically by 0.5 for 
display purposes. 



2/3 level, by a-elements. We also present an average spec- 
trum of ellipticals/spheroids, built from a homogenized sub- 
set of our sample that can be used as template of the inte- 
grated early-type galaxy spectra. 
4. The strength of Mgl feature at 1 .5 yum was measured for the 
first time in a heterogeneous sample of galaxies, providing a 
new constraint for future NIR spectral synthesis models. 
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